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ABSTRACT 

Weak lensing by large-scale structure provides a direct measurement of matter fluctuations in the 
universe. We report a measurement of this 'cosmic shear' based on 271 WFPC2 archival images from 
the Hubble Space Telescope Medium Deep Survey (MDS). Our measurement method and treatment of 
systematic effects were discussed in an earlier paper. We measure the shear variance on scales ranging 
from 0.7' to 1.4', with a detection significance greater than 3.8cr. This allows us to measure the nor- 
malization of the matter power spectrum to be as = (0.94 ± 0.10 ± 0.14)(0.3/^^m)°■'"'(0.21/^)°•l^ in a 
ACDM universe. The first la error includes statistical errors only, while the latter also includes (gaus- 
sian) cosmic variance and the uncertainty in the galaxy redshift distribution. Our results are consistent 
with earlier cosmic shear measurements from the ground and from space. We compare our cosmic shear 
results and those from other groups to the normalization from cluster abundance and galaxy surveys. 
We find that the combination of four recent cosmic shear measurements are somewhat inconsistent with 
the recent normalization using these methods, and discuss possible explanations for the discrepancy. 

Subject headings: cosmology: observations — dark matter — gravitational lensing — large-scale structure 
of the universe 



1. INTRODUCTION 

Weak gravitational lensing by large-scale structure has 
been shown to be a valuable method of measuring mass 
fluctuations in the universe (see Mellier at al. 2001 for 
a review). This effect has been detected both from the 
ground (Wittman et al. 2000; van Waerbeke et al. 2000, 
2001; Bacon et al. 2000, 2002; Kaiser et al. 2000; Hoek- 
stra et al. 2002) and from space (Rhodes, Refregier, & 
Groth 2001, RRGII; Hammerle et al. 2001). These results 
bode well for the prospect of measuring cosmological pa- 
rameters and the mass distribution of the universe using 
weak lensing. 

In this letter, we present the highest significance de- 
tection of cosmic shear using space-based images. It is 
based on images from the Hubble Space Telescope (HST) 
Medium Deep Survey (MDS; Ratnatunga et al. 1999). 
We apply the methods for the correction of systematic 
effects and detection of shear we have previously devel- 
oped (Rhodes, Refregier, and Groth 2000; RRGI) to 271 
WFPC2 fields in the MDS. The method of RRGI is specif- 
ically adapted to HST images and takes advantage of the 
small PSF of the HST (0.1" as compared to ~0.8" from 
the ground) . This affords us a higher surface density of re- 
solved galaxies as well as a diminished sensitivity to PSF 
smearing when compared to ground-based measurements. 
We develop an optimal depth- weighted average of selected 
MDS fields to extract a weak lensing signal. We then use 
this signal to derive constraints on the amplitude of the 
mass power spectrum and compare this to measurements 
from previous cosmic shear surveys and from other meth- 
ods. 

2. DATA 
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The MDS consists of primary and parallel observations 
taken with the Wide Field Planetary Camera 2 (WFPC2) 
on HST. We selected only the I-band images in chips 2,3, 
and 4 to study weak lensing. To ensure random lines-of- 
sight, we discarded fields which were pointed at galaxy 
clusters, leaving us with 468 I-band fields. To ensure that 
our fields are independent, we selected 291 fields separated 
by at least 10', beyond which scale the lensing signal drops 
considerably (see Figure 2). 

We used the MDS object catalogs (Ratnatunga et al. 
1999) to determine the position, magnitude, and area of 
each object, as well as to separate galaxies from stars. 
We used the chip-specific backgrounds listed in the MDS 
skysig files, which are consistent with backgrounds cal- 
culated using the IRAF task imarith. Not using object- 
specific backgrounds necessitated the discarding of another 
20 fields with a large sky gradient. Our final catalog thus 
consisted of 271 WFPC2 fields amounting to an area of 
about 0.36 deg^. 



3. PROCEDURE 

The procedure we used for measuring galaxy ellipticities 
and shear from the source images is described in detail in 
RRGI (1999) (see also RRGII and Rhodes 1999). It is 
based on the method introduced by Kaiser, Squires, and 
Broadhurst (1995), but modified and tested for applica- 
tions to HST images. The usefulness of our method was 
demonstrated by our detection of cosmic shear in the HST 
Groth Strip (RRGII). 

We correct for camera distortion and convolution by the 
anisotropic PSF using gaussian-weighted moments. Cam- 
era distortions were corrected using a map derived from 
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stellar astrometric shifts (Holtzman, et al., 1995). PSF 
corrections were determined from HST observations of four 
stellar fields These fields were chosen to span the focus 
range of the HST as shown by Biretta et al. (2000). Fi- 
nally, we derive the ellipticities e,; of the galaxies and con- 
vert them into shear estimates using 7^ — G^^ei, where G 
is the shear susceptibility factor given by equation (30) in 
RRGI. 
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Fig. 1. — Distribution of the median magnitudes Im, of the fields 
(soHd line and left j/-axis). The top x-axis shows the approxi- 
mate corresponding median redshift according to equation Ivt). The 
weights of the fields (ui j ) averaged in each magnitude bin are shown 
as the dashed line and right j/-axis, in arbitrary units. The eff'ectivc 
depth of the survey with this weighted scheme is shown as the ver- 
tical dotted line and corresponds to Im = 23.5 or ZfYi = 0.95 ±0.10. 

To limit the impact of noise and systematics, we made 
a number of cuts to select our galaxy sample. We first 
discarded objects for which e > 4 after deconvolution. To 
avoid using galaxies with low signal-to-noise, we included 
only galaxies that have a magnitude / < (J^ -I- 2), where 
is the median magnitude (before the cut) of the field 
and chip in which the galaxy lies. This is consistent with 
the magnitude cut we made in RRGII. We also discarded 
small galaxies (size d < 1.5 pixels) in order to minimize the 
effects of the anisotropic PSF on our measurements (see 
RRGI). The final galaxy sample contained about 3.1 x 10^ 
galaxies. The distribution of the median magnitude Im 
(after the cuts) of our fields is shown in Figure 1, which 
reveals a wide range of depth for the MDS. 

For the magnitude range of the MDS (19 < / < 27), wc 
use spectroscopic redshift determinations from the DEEP 
survey (DEEP Collaboration, 1999) and the Hubble Deep 
Fields (Lanzetta et al. 1996) to determine that the median 
redshift Zm in a field is related to the median magnitude 
/,„ by the equation 



0.722 + 0.149 X (/„ -22.0) 



(1) 



This agrees well with an extrapolation of the CFRS red- 
shift distribution (Lilly et al. 1995) we used to determine 
the median redshift of objects in the Groth Strip (RRGII). 
Both methods give Zm = 0.9 ± 0.1 for the Groth Strip 



(/„j = 23.6), where the error gives a measure of the la 
systematic uncertainty in the above relation. 



4. ESTIMATOR FOR THE SHEAR VARIANCE 

We wish to derive a measure of the shear variance on 
different angular scales by averaging over the Nc — 3 
chips in each of the Nf = 271 fields. As explained in 
the fields have varying depths but are are sufficiently 
far apart to be statistically independent. As in RRGII, 
the total mean shear 7ic/ in chip c and field / can be 
measured by averaging over all the selected galaxies that 
it contains. It is equal to the sum of contributions from 
lensing, from noise and from systematics, and can thus be 
written as jicf — jj^Y + licf^'^ + 7iJ/- ^he noise vari- 
ance cTnoisc c/ = {hcT^°\'^) measured from the data 
by computing the error in the mean 7ic/ from the dis- 
tribution of the galaxy shears inside the chip. As RRGI 
showed, the systematics are greatly reduced if the shear 
is averaged over the chip scale and if small galaxies (with 
d < 0.15") are discarded. In this case, the systematics are 
dominated by the time-variations of the PSF and induce 
a shear variance a'^^^ 
0.0011^ (see RRGII). 



(l7cFI^) approximately equal to 



For each field /, an estimator for the shear variance uj^j^g 
on the chip scale is given by 
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where cr^^i.e,/ ^N^^J^c 



c noise, c/ ' 



and similarly for Cgyg ^ . 



Assuming gaussian statistics, the error variance of the 
combined estimator is given by (see Bacon, Refregier & 
EUis 2000 for the case = 2) 
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where CT^^^,^ ^ ^^lens^lenp ^ ^^lens^lens^ jg ^j^g 

correlation between chips c and c' and can also be mea- 
sured from the data. The term crj^^ns f corresponds to cos- 
mic variance, and the last term arises because the chips 
within a field are not statistically independent. While the 
lensing shear field is known to be non-gaussian on scales 
smaller than about 10' (eg. Jain & Seljak 1997), the non- 
gaussian corrections to this error estimate are small for 
noise-dominated 2-point statistics like the variance (see 
discussion in RRGII and White & Hu 2000). 

Because the fields have a range of depths, it is desir- 
able to combine the individual estimators cr^iens,/ using a 
weighting scheme of the form 



Cr^lens = Wf aliens, f 




(4) 



A convenient choice for the weights is given by 

= <o1se,f> (5) 
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i.e. to the inverse square of the noise contribution to the 

error in cr^ions,/- This weighting scheme is nearly optimal, 
and avoids including the lensing signal crf^^-^^ j itself. The 
average weights {wf) in several magnitude bins are shown 
in Figure 1 . As expected, deeper fields have larger weights 
since they contain a larger number of galaxies and thus 
have a smaller value of (Tnoisc,f • 

To measure the shear variance on the field scale, we first 
average the shear within each field and apply the same pro- 
cedure. This time, however, the cross-correlation term in 
Equation vanishes since each field is independent. Sim- 
ilarly, we can consider pairs of chips to measure the shear 
variance on intermediate scales. 
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Fig. 2 

circular celh Our observed value (HST) as well as that observed by 
other groups: van Waerbeke et al. (2001, CFHT vW+), Wittman 
et al. (2000, CTIO), Kaiser, Wilson & Luppino (2000, CFHT K-I-), 
Maoli et al. (2000, VLT), Bacon et al. (2002, WHT and Keck). For 
our measurement, the inner error bars correspond to noise only, while 
the outer error bars correspond to the total error (noise + cosmic 
variance). The errors for the measurements of Maoli et al. (2000) 
and van Waerbeke et al. (2001) do not include cosmic variance. The 
measurements of Hammerle et al. (2001) and Hoekstra et al. (2002) 
are not displayed but are consistent with the other measurements. 
Also displayed are the predictions for a ACDM model with Q,m = 0.3, 
erg = 1, and F = 0.21. The galaxy median redshift was taken to be 
= 1.0,0.9, and 0.8, from top to bottom, respectively. 



this Eq.). For instance, on the chip scale we obtain 
o'iens(0-''2') = (3.5 ± 0.9 ± l.f) X fO-^ yielding a de- 
tection significance (inner error) of 3.8tT with this scale 
alone. As a check of systematics, we analyzed our signal 
into E and _B-modes using the aperture mass Map (0.67') 
statistic on the chip scale (Schneider et al. f998; van 
Waerbeke et al. 2001). For i<^-modes, we find the up- 
per limit {Ml^ j^) = (0.4 ± f.7) X 10-4 (Icr) which is 
consistent with the signal expected in a ACDM model 
((Mfp_^) ~ 0.6 X lO"'^; Schneider et al. 1998). For the 
B-modes, we find {Ml^^g) = (0.3 ± 1.7) x lO"'^ (Id), as 
expected in the absence of systematics (which corresponds 

to (A/fp^^) EE 0). 

The measurements from other groups are also plotted 
in Figure 2, along with the predictions for a ACDM model 
with (Tg = 1 and F — 0.21. The central value for F is 
close to the recent measurement of this parameter from 
galaxy clustering (eg. Percival et al. 2001), while keeping 
the F = Vimh relation valid for h — 0.7. The predictions 
are plotted for a range of galaxy redshifts Zm = 0.9 ±0.1, 
corresponding approximately to the uncertainty and dis- 
persion of this parameter in the different surveys. In our 
case, the effective median I-magnitude of our measurement 
is Im = '^fWflmj/^fWf ~ 23.5, which corresponds to 
a median redshift of = 0.95 ±0.10 (see Eq. ||l|). The ef- 
fective magnitude and corresponding redshift are plotted 
in Figure 1. Given the range of median redshifts in the 
different surveys and the correlation between angular bins 
for the variance, our results are in good agreement with 
these other measurements and with the ACDM model. 

Our measurements can be used to constrain cosmolog- 
ical parameters. Because our measurements on different 
scales are not independent, we conservatively only consider 
the shear variance on the chip scale [9 = 0.72'). Within 
a ACDM model, it is predicted to be (see RRGII), within 
an excellent approximation. 
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where erg is the amplitude of mass fluctuations on 
8 ft,-^Mpc scales, and Vtm is the matter density pa- 
rameter. Inverting this equation, we find that our 
measurement of (yfc^g yields CTg = (0.94 ± 0.10 ± 
0.12)(r!„/0.3)-"-44(r/0.21)-0-i5(^„/0.95)-o™, where 
the first error is statistical only and the second also in- 
cludes cosmic variance. To this error must be added that 
arising from the uncertainty in the median effective red- 
shift Zm = 0.95 ± 0.10. After propagating this error, we 
obtain 



5. RESULTS 

Our measurement for the shear variance cfi^^s{9) on dif- 
ferent scales is shown in Figure 2. The angular scale 
is the radius of an effective circular cell whose mean pair 
separation equals that of the chip configuration consid- 
ered {9 ~ 0.72', 1.11', and 1.38', for 1, 2 and 3 1.27' chips, 
respectively). The outer Icr error bars include both sta- 
tistical errors and cosmic variance (from Eq. |^]), while 
the inner error bars only include statistical errors (i.e. 
by setting rjy^^^ and ct^ to on the right-hand side of 



CT8 = (0.94 ±0.10 ±0.14) ( ^ 



0.44 



0.21 



0.15 



(7) 



where the first error reflects statistical errors only, and 
the latter is the total error and includes statistical errors, 
cosmic variance, and redshift uncertainty (all Ict). 

Figure 3 shows the comparison of our measurement of trg 
(IIST/WFC2) with that from other weak lensing surveys 
and from other methods. A ACDM model with $7™ — 0.3 
and F = 0.21 was assumed (except for van Waerbeke et al. 
2001 who marginalized over F). Our erg value is consistent 
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with the other recent cosmic shear measurements of Bacon 
et al. (2002), Hoekstra et al. (2002), and van Waerbeke 
et al. (2002) and also with the 'old' normalization from 
cluster abundance (eg. Pierpaoli et al. 2001). This was 
recently revised to a lower normalization, by using the ob- 
served mass-temperature relation rather than that derived 
from numerical simulations (eg. Seljak 2001). A similar 
normalization was derived by combining measurements of 
galaxy clustering from 2dF and CMB anisotropy (Lahav 
et al. 2001). Our results are consistent with this new 
normalization at the 1.4cr level. 

ACDM, n„=o.3, r=o.2i 
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Fig. 3. — Comparison of the determination of trg by different groups 
and methods. The errors have all been converted to Icr, and a ACDM 
model with Qrn = 0.3 and F = 0.21 was assumed (except for van 
Waerbeke et al. 2002 who marginalized over F between 0.1 and 
0.4). The vertical dotted lines shows the weighted average (weights 
oc a-~^) of the 4 cosmic shear measurements and associated la error. 



6. DISCUSSION AND CONCLUSION 



We have achieved the highest significance detection of 
cosmic shear using space-based images to date. Using the 
MDS, we have detected the shear variance on 0.7' to 1.4' 
scales with a significance greater than 3.8(J. From the 
amplitude of the signal we derived a normalization of the 
matter power spectrum given by Equation ^ , with errors 
which include statistical errors, (gaussian) cosmic variance 
and the uncertainty in the galaxy redshift distribution. 
Our results agree with previous measurements of the rms 
shear using both ground and space based images at the 
Icr level and with the 'old' (eg. Pierpaoli et al. 2001) and 
'new' (eg. Seljak 2001) cluster abundance normalization 
at the 0.4cr and lAa level, respectively. 

A weighted average of the four recent cosmic shear mea- 
surements shown in Figure 3 yields erg — 0.91 ± 0.02, for 
a ACDM model with n„i = 0.3 and F = 0.21 (see ver- 
tical bars in Figure 3). This is consistent with the old 
cluster normalization at the 1.2cr level, but somewhat in- 
consistent with the new cluster normalization at the 2.5a 
level, where the uncertainty is dominated by that from 
cluster abundance. This discrepancy could be caused by 
unknown systematics in the cluster abundance or cosmic 
shear methods. For the latter case, the calibration of the 
shear measurement methods would need to be revisited, in 
the context of current and upcoming surveys. The inaccu- 
racy of the calculation of the non-linear power spectrum 
and of the halo mass function may also contribute to the 
error budget. If confirmed however, this discrepancy could 
have important consequences for our understanding of the 
physics of clusters, or may require extensions of the stan- 
dard ACDM paradigm for structure formation. 

We thank David Bacon, Richard Massey and Richard 
Ellis for fruitful discussions, and the anonymous referee 
for useful suggestions. AR was supported by an EEC fel- 
lowship from the TMR network on Gravitational Lensing 
and by a Wolfson College Research Fellowship. EJG was 
supported by NASA Grant NAG5-6279. JR was supported 
by an National Research Council-GSFC Research Associ- 
ateship. The Medium Deep Survey catalog is based on ob- 
servations with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research 
in Astronomy. 



REFERENCES 



Bacon, D.J., Refregier, A., & Ellis, R.S, 2000, MNRAS, 318, 625 
Bacon, D.,J., Massey, R., Refregier, A., & Ellis, R.,2002, submitted 
to MNRAS 

Biretta, J. A., et al. 2000, WFPC2 Instrument Handbook, Version 

5.0 (Baltimore: STScI). 
DEEP Collaboration, 1999, Private Communication 

Hammerle et al., 2001, A&A, 385,743 . . 

Hoekstra, H. et al., 2002, ApJ in press, preprint astro-ph/0202285 

Holtzman, J. et al., 1995, PASP, 107, 1065 

Jain B., Seljak U., 1997, ApJ, 484, 560 

Kaiser, N., & Squires, G., & Broadhurst, T. 1995, ApJ, 449, 460 
Kaiser, N., Wilson, C, & Luppino, G.A., 2000, submitted to APJL, 

preprint astro-ph/000338 I . 

Lahav, C, et al. 2001, submitted to MNRAS, Eistro-ph/0112162 
Lanzetta, K., Yahil, A , ie Frrriandn7-Sntn, A , I MMbT4Nature7a81, 



759 and data at 
Lilly, S., et al. 1995' 
Maoli, R., et al. 2001 
Mellier, Y., et al.. 



ahiJ 
,t \i 

995 



ttp: //www. ess . sunysb . edu/astro/ 
TtpT 



4bb, 108 
A&A, 368, 766 
2001, Deep Fields 
ESO/ECF/STSCI Workshop, Springer 



Proceedings of the 



Percival, W., et al. 2001, MNRAS, 327, 1297 
Pierpaoli, E., Scott, D., & White, M., 2001, MNRAS, 325, 77 
Ratnatunga, K., Griffiths, R., & Ostrander, E., 1999, AJ, 188, 86 
Rhodes, J. 1999, PhD Thesis, Department of Physics, Princeton 
University 

" ~ " & Groth, E.J., 2000, (RRGI), ApJ, 536, 



Groth, E.J., 2001 (RRGII), ApJ, 552L, 



Jain, B., & Kruse, G., 1998, MN- 



Rhodes, J., Refregier, A 
79 

Rhodes, J., Refregier, A. 
85 

Schneider, P., van Waerbeke, L 

RAS, 296, 873 , . 

Seljak, U., 2001, submitted to MNRAS, preprint astro-ph/011136; 

van Waerbeke et al., 2000, A&A, 358, 30 

van Waerbeke et al., 2001, A&A, 374,757 . 

^^^^n Waerbeke et al., 2002, submitted to A&A, preprint astro- 



ph/0202503 
White, M., & Hu, W., 2000, ApJ, 537, 1 

Wittman, D.M., Tyson, J. A., Kirkman, D., Dell'Antonio, I. 
stein, G. 2000, Nature, 405, 143 



Bern- 



